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Abstract 

We investigate correlations among various parameters in the solution space of mini- 
mal supersymmetric grand unification. In particular the extent to which the top quark 
Yukawa coupling exhibits fixed point behavior is discussed and we compare various an- 
alytic approximations to its value at the top mass with its exact value in numerical 
solutions. 

One of the successes of the idea of supersymmetric grand unification is the prediction of the b quark to 
r lepton mass ratio and the accompanying prediction of a top quark mass significantly above that of the Z. 
These predictions are based on the solution of a set of coupled renormalization group differential equations 
involving the gauge and Yukawa couplings. At least in the case of small tan(/3), (< 5), the solution space 
of minimal SUSY unification is a ten dimensional space defined by the values of the following ten parameters. 

1) A unification scale Mx 

2) A unified gauge coupling ao(Mx) 

3) A top Yukawa at M x , a t {M x ) 

* Talk presented at the Workshop on Yukawa Couplings and the Origin of Mass, Gainesville, Florida, Feb. 
11-13, 1994 



4) A Susy scale M s 

5) A ratio tan(/3) of the Higgs vacuum expectation values 

6) The weak angle sin 2 (9w) 



7) The fine structure constant a(Mz) 

8) The strong coupling constant a^{Mz) 

9) The value of the top quark mass M t 

10) The value of the b/r mass ratio at the b quark scale 

A " solution" is defined as a set of values for these ten parameters which is consistent with the renormal- 
ization group running and with the experimental constraints: 

a(M z ) = 127.9 ± 0.2 (1) 
sm 2 {0 w {M z )) = .2328 ± .0007 (2) 

m b /m T = 2.39 ±0.10 (3) 

The latter corresponds to a physical b quark mass of 4.95 ± .15GeV. There are some who feel that the 
uncertainty in this quantity is much smaller than taken here but out of respect for the complications of 
confinement we content ourselves with this 4.5% uncertainty. The uncertainties on the other two quantities 
are below 1% and once the top quark mass is known the current data will specify them to about 0.15% due 
to the correlation 

sm 2 (6 w ) = .2324 - .002(M 2 /(138Gey) 2 - 1) ± .0003 (4) 

The hope is that the increasing precision with which these numbers are known will shed light on the physics 
at the scales M$ and Mx and will enable predictions to be made for mt and a^{Mz)- 

In addition to the experimental constraints eqs. 0J|^ we assume that ctt(Q) < 1 for all Q ("perturbativ- 
ity") and that 

WOGeV < M s < M s , max (5) 

If the Susy scale were below lOOGeV, with the expected degeneracy splittings among the Susy particles, we 
would have expected unacceptably large contributions to the Z width due to Susy decay modes. If Ms is 
too large the theoretical benefits of Susy in explaining the stability of a low scale for electroweak symmetry 
breaking is lost. Ms,max is variously taken to be ITeV or lOTeV. Ideally this argument would prefer Ms 
below one TeV since the electroweak scale is in the hundred GeV range. In addition, a Susy solution of the 
dark matter problem would require the mass of the LSP to be no higher than 200Gey again suggesting an 
average Susy mass below one TeV. Nevertheless we will, for the sake of conservatism, take Ms, m ax — lOTeV. 
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In the light gluino scenario, tan(/3) is restricted to be between 1 and 2.3 m and, for simplicity, we limit 
our investigation in the heavy gluino case, to the range 1 < tan(/3) < 5 so that we can neglect the effect 
of the b Yukawa on the running of the couplings. This range is also preferred by proton decay. Radiative 
electroweak breaking would predict a value of tan(/3) very close to 1.8 in the light gluino case. 

In the simplest version of SUSY unification one assumes that all the GUT scale particles are degenerate 
at Ms and all the SUSY partners of the standard model particles are degenerate at Ms- In the light 
gluino variant one assumes that the partners of the squarks and sleptons are at Ms together with the heavy 
Higgs, while the photino and gluino are in the low energy region below Mz and the other neutralinos and 
charginos are at the scale of Mz as suggested by the Mi/ 2 = model. This is consistent with the current 
experimental gluino searches which leave open (at least) the three windows shown in fig.l. Fig.l updates 
the chart published by the UA1 M group in 1987 to include the LEP results which are probably the most 
model independent constraints together with the results of the HELIOS || collaboration which searched 
for weakly interacting neutral particles. 

In higher level variations, the GUT scale spectrum is assumed to be non-degenerate and possibly richer 
than in the minimal model and/or the Susy scale Ms is split into different masses for the various particles. 
In the latter schemes for each non-degenerate spectrum of Susy particles there is an effective degenerate 
scale Ms which leads to the same unification solution apart from small two loop effects. Q Each solution in 
the degenerate case corresponds to a family of solutions with different splittings among squarks and sleptons 
and, as long as Ms is above Mz , no solutions for the ten parameters above are lost by assuming degeneracy. 

The two-loop differential equations to be solved are summarized in the papers of re/s.|^, |(|. We follow 
a "top-down" approach where one begins by choosing random values for the first 5 quantities in the list 
above, extrapolating to low energies and discarding the choices which are inconsistent with the constraints 
of egsj^^,^, ||. The surviving solutions are stored in a data set that can be queried for correlations among 
the various parameters. The solution space forms a small connected region in the ten dimensional space. 
Careful checking is required around the borders of this region to insure that these are in fact the borders and 
that no solutions exist outside this region. In fig.2a for example we show the solution space of the minimal 
Susy model projected onto the Ms — Mx plane in a random non-overlapping sample from 813 solutions 
while fig.2b shows the same projection for a non-overlapping subset of 1200 solutions in the light gluino 
scenario. The shape coding labels the a${Mz) value of each solution. In the heavy gluino case solutions are 
found only for 0.111 < a^(Mz) < -134 while in the light gluino case the range is 0.122 < a^{Mz) < .133. 
These ranges are divided into quadrants indicated in the solutions of fig.2 by rectangles, triangles, ovals, 
and diamonds for a.${Mz) in the lowest to the highest quadrant respectively. One sees that with the light 
gluino option, The GUT scale is restricted to values comfortably above 10 16 GeU while this is not the case 
in the standard Susy picture. In Susy unification, proton decay via lepto-quark gauge bosons of mass as low 
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as 10 15 GeV is not in contradiction with current limits. However, proton decay via the super-heavy scalars 
requires these particles to have masses in excess of 10 16 GeV. Since the GUT scale, above which the theory 
is grand-unified, is the maximum mass of the GUT scale particles, solutions with Mx below 10 16 GeV are 
probably not acceptable. One could therefore add a sixth experimental constraint 

M x > 10 w GeV (6) 

which would cut the solution spaces of figs.2a, b at the corresponding limit. The solution sets of fig.2a, b 
were generated in the course of work reported in ref. [Q. One sees in figs.2a, b the effect reported there that 
the «3 values in minimal Susy unification are significantly more constrained in the light gluino scenario than 
in the usual picture. This result, however, is critically dependent on the super-gravity inspired prediction that 
the charginos and neutralinos are relatively light (below the Z) when the gluino is light. In these data sets, 
the constraint of eg.|]was not enforced, although that of eg.||was, so that a Susy prediction for sin 2 (#vK) could 
be made. This prediction can be read from figs.3a,b which show the solutions in the sin 2 (#vv) — a^{Mz) 
plane. One sees that in both heavy and light gluino cases, the predicted values of sin 2 (9w) he between 
0.230 and 0.2325. Thus the minimal model predicts the weak angle with a 1% accuracy and agrees with 
experiment. The shape coding indicates the quadrant values of Ms in the range from lOOGeU to lOTeU 
(rectangle, triangle, oval, diamond from lowest to highest quadrant). In figures 4a, b we show the heavy and 
light gluino solutions respectively in the tan(/3) — M t plane with the values of Ms indicated by the shape 
coding. One sees that in both cases the top quark mass is bounded below by about 143GeU and that in the 
light gluino case the band of solutions is appreciably broader. It is interesting to note that the preliminary 
evidence from Fermilab for M t ~ HAGeV suggests a tan(/3) ~ 1.8 as required in the light gluino scenario 
with radiative electroweak breaking. Q In the current work, however, and that of |Q, radiative breaking is 
not assumed. 

Much has been written |^, |l(| about the quasi-fixed point behavior of the top Yukawa. This can be 
roughly defined by the statement that the value of the top Yukawa at the top mass given by 

a t (M t ) = Mf • (173GeUsin(/3)(l + 4a 3 (M t )/(37r) + ll(a 3 (Af z )/^) 2 ))- 2 (7) 

is dependent only on the gauge couplings at the top scale and is independent of the GUT scale parameters at 
least for some range of those parameters that is consistent with the Susy unification solutions. Our present 
purpose is to clarify and quantify this statement. 

Neglecting the effect of the b and r Yukawas on the running of a t and two loop contributions, the top 
Yukawa in the Susy region satisfies 

2irda t /dt = a t {Qa t - ci^oti) (8) 
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where the at are the three gauge couplings, t = ln(Q), and 

c M = (13/15,3,16/3) (9) 

for i = 1,2,3 in that order. If at is initially (at the GUT scale) higher than Ct^ai/6 the naive prediction is 
that it will fall with decreasing Q until the right hand side of eg.|^ vanishes. If at is below ct.iO-i/Q initially, 
it will rise toward that value as Q decreases. At the most naive level, one might expect from eqJ^ that at(Q) 
will approach the fixed point expression 

a ( 3(Q)=c t>i a i /6 (10) 

From our unification solution set it is a simple matter to calculate for each solution the top Yukawa at 
M t from eg.J^. The gauge couplings at M t can be related to the ai{Mz) from the renormalization group 
expressions. We can therefore check how well at approaches egjl^ at Q = M t . In figures 5a, 6, for the 
heavy and light gluino scenarios respectively, we show the correlation between at(M t ) and a^j{Mt). The 
correlation is far from the close equality one might have expected. The shape coding, rectangles, triangles, 
ovals, diamonds, indicates an Ms value in the lowest to the highest quadrant respectively with the total 
range being lOOGeV < M s < lOTeV. We will return later to discuss the tail of solutions out to high 
a^f(M t ) evident in figures 5a, 6. 

Ibanez and LopezQ have shown that eg. 8 is satisfied by a top Yukawa given by 

-a t {M x )F{Q) 

at{Q) =l + 6at(M x )F(Q)/4n (U) 

where 

F(e) = -£f -rf-jCs^W" (12 > 

and 

HQ) = Q^q f ^Q) ( 13 ) 

.Eg. 11 is of course not an analytic solution of eg.|] since the Q' integral in eg. 12 is not analytically soluble. 
However it does illustrate the quasi- fixed-point behavior in that at(Q) becomes independent of at{Mx) if 
F(Q) is sufficiently large as would happen, for example, for large enough Mx- From eg.|ll] one can write a 
second expression for the fixed point by dropping the 1 in the denominator of egjll|. 

^ = ^4 lnF(0) (14) 

Because we don't have here an analytic solution for F(Q), it is difficult to predict the extent of the inde- 
pendence of a t (M t ) on Mx and a t (Mx)- Eqn4 does indicate dependence of a t (M t ) on Mx and indirectly 
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on at(Mx) through the two loop effects on the running of the gauge couplings. The empirical dependence 
of at(Mt) on Mx is substantial as can be shown by projecting the solution space onto the ctt(Mi) — Mx 
plane. The Ibanez-Lopez expression for the quasi-fixed-point does not yield analytically the dependence of 
at(M t ) on the ten parameters of the unification solution. We seek preferably a quasi-fixed-point expression 
analogous to that of eg.|lO|. 

One of the problems with both of the above treatments is that the running of the top Yukawa changes 
dramatically below the Susy scale. Then instead of eqs. | and | we have 

dry sm 9 

27T —jj- = a t \~2 a t ~ c t,i a *) ( 15 ) 

where below M$ one defines the effective top Yukawa 

a s t m = a t sin 2 (16) 

and 

c s t ™ = (17/20,9/4,8) (17) 



If Ms is above M t one extrapolates from Ms to M t using egs.|l5| and 17 and then redefines at(M t ) through 
eg.|l6| before substituting in eg.f^. This rescaling has only a higher order effect on the final answer so we use a 
rescaling by the fixed value sin/3 = .7641 in both the numerical running and the analytic expressions below. 

Thus, if the top Yukawa is following a quasi-fixed-point given approximately by eg.[l0]down to Ms, below 
the Susy scale we would expect it to be drawn toward a naive fixed point corresponding to 

°%f m (Q) = \<7 ( 18 ) 



The corresponding otpj is more than twice the fixed point of the Susy regime suggested by eg. 10. This 



effect could partially explain the large deviations of at(M t ) from a^/AMt) evident in figures 5a, 6. However, 
it is clear that even with Ms as large as WTeV, the top Yukawa does not have time to reach a standard 
model quasi- fixed-point. From figures 5a, b one sees a tendency for the top Yukawa to rise with increasing Ms 
as would be expected from this effect but it never achieves the doubling expected naively from egjl^. Thus, 
if Ms is in the 1 to lOTet^ region, the top quark Yukawa is unlikely to have reached a limiting behavior. The 
tail of events at large ce^j in figure 5 represents the events in which Ms < M t so the Susy quasi-fixed-point 
should be most accurately attained. Paradoxically, it is here that the discrepancy between at(M t ) and a^J 
is largest. Clearly a more accurate representation of the top Yukawa behavior is required. 
In the SUSY regime the gauge couplings change according to the law 

2^ = a% (19) 
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with 

bt = (33/5, 1, -3) (20) 

We can combine eqs.^ and |l9| with arbitrary parameters Xi to write 

d(a t - Ai«i) 2 2 
27T — = &a t - ct^otiOtt - Aia^i (21) 

where summation over repeated indices is intended and, for the purpose of an analytic approximation, we 
have again neglected two loop contributions to the running. Ideally, one would seek such that when at 
reaches 

a t; f = Xiai (22) 

the right hand side of eq.^l] would vanish. Thus we would seek solutions to the equation 

6(Aiai) 2 - i-!,,n,X,n , - Xia?bi = (23) 

However, no set of exists that will satisfy eg.^3] for arbitrary values of the c^. We propose therefore to 
write eg.^l] in the form 

ir d(a t -Xtai) 



dt 



where 



We will now choose 



(a t - ct^/12) 2 - S 2 (24) 



(c Mai /12) 2 + XAa 2 /6 (25) 



Xi = c M /12 (26) 



and make the approximation of ignoring the Q dependence of 8. Then eg|2J can be integrated to write 

at (Q) = Kai (Q) + % + (5 _(, _^MQ/M,)/, ( 27 ) 

where yo = at(Mx) — Xiai(Mx) and 6 is defined at Q. This is an identity at Q — Mx and gives an excellent 
approximation to at(Q) for all Q > Ms ■ In the limit Mx/Q — > oo, the exponentials in eg.|2^ become 
negligible and at(Q) is determined solely by the ai(Q) becoming independent of the GUT scale values. In 
this case one could talk about a quasi-fixed-point behavior. However in practice Mx is never large enough to 
justify dropping the exponentials. Furthermore, when M t < Ms we must face the complication of integrating 
at in the standard model region. In this case we use et?.^ for a t (Mg) and write the approximate form using 

a t (M t ) = at(M s ) (l + lnM ^ Ms ( ,7641 • % -a t {M s ) - dgai(M t jf) (28) 
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(3) 

In /zg.6a, b we show the correlation between the right hand side of this equation which we may call a t .^ 
and the exact numerically calculated left hand side. The prediction is seen to hold within a few percent 
over the entire solution space. Given M t and Ms , the gauge couplings at those scales can be found to 
sufficient accuracy by the first order extrapolation from their values at Mz ■ If the exponentials in e<?.p7| 
were negligible, egs.f?] and 28, together with knowledge of M t , A/5 and the gauge couplings at the Z would 
yield an accurate measure of tan(/3) independent of the GUT scale parameters. However, only to a crude 
approximation ~ 20% can the exponentials in eqX2l\ be neglected. 



The author is indebted to Philip Coulter for discussions in the course of this work. The research reported 
here was supported in part by the Department of Energy under Grant No. DE — FG05 — 84i?i?40141. 
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FIGURE CAPTIONS 

Fig.l. Low mass windows for the gluino mass as a function of the squark masses. The hatched 
areas are disfavored by the indicated experiments. The dot-dashed curves represent the loci of ex- 
pected gluino lifetimes 10~ 6 ,s, lCT 8 ,s, 10 _10 s, and 10 _12 s respectively from the highest to the lowest 
curve. 

Fig.2. (a) The correlation between the SUSY scale, Ms, and the GUT scale, Mx, in the heavy 
gluino case, m g = Ms- Each solution corresponds to an allowed point in the ten dimensional space 
discussed in the introduction. For given Ms, no solutions exist outside of the broad band shown. 
The width of the band corresponds to summing over all other eight parameters. The aj,{Mz) value 
for each solution is indicated by the shape coding (see text). 
(b) The same correlation in the case of the light gluino (nig < Ms). 

Fig.3. (a) The correlation between a>t(M t ) and sin 2 (9w) m the heavy gluino case. Solutions in 
the 1st, 2nd, 3rd, and 4th quadrant of the Ms range, lOOGeU < Ms < lOTeU, are printed as 
rectangles, triangles, ovals, and diamonds respectively. 
(b) The same as (a) but for the light gluino scenario. 

FigA. (a) The correlation between tan(/3) and M t in the heavy gluino case. The Ms quadrants 
are indicated as in Fig.3. (b) The same as (a) but for light gluinos. 

Fig.5. The correlation between the "naive" top Yukawa fixed point a^j(M t ) and the actual top 
Yukawa a>t(M t ) found in the numerical solutions in the heavy gluino case (a) and in the light gluino 
case (b). 10% to 20% departures are observed. The quadrant values of Ms are indicated by the 
shape coding as in fig.3. 

(3) 

Fig. 6. The correlation between the approximate analytic value, a^.j(Mf), of the top Yukawa and 
the exact numerical value of at(Mt) in the heavy gluino solution space (a) and the light gluino 
solution space (b). Agreement within 2% is found. The quadrant values of Ms are indicated by 
the shape coding as in fig.3. 
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